Surface-induced layer formation in polyelectrolytes 
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We analyze, by means of an RPA calculation, the conditions under which a mixture of oppositely 
charged polyelectrolytes can micro-segregate in the neighborhood of a charged surface creating a 
layered structure. A number of stable layers can be formed if the surface is sufficiently strongly 
charged even at temperatures at which the bulk of the mixture is homogeneous. 



I. INTRODUCTION 

It has been described in an article by Decher ||,|2j a 
procedure for creating a coating for a (charged) surface 
that consist of alternating layers of two different, oppo- 
sitely charged, polyelectrolyte homopolymers. The prop- 
erties and uses of these systems has been explored by 
several groups 0-0]. The basic technology for the cre- 
ation of these coatings is easy to describe. The starting 
surface is charged (say, positively) and it is first dipped 
into a solution of negatively charged homopolymers and 
then withdrawn from it. This forms a first layer since 
the polyelectrolyte is strongly attracted to the surface. 
The surface is washed and one gets rid of excess polymer 
chains that are only loosely tied to the surface. The sur- 
face is now dipped in a solution of the second homopoly- 
mer that is positively charged. This second homopolymer 
is strongly attracted to the coated surface, which is now 
effectively negatively charged. The process is repeated 
several times, and the layering is found to be stable. A 
more detailed exposition is presented in the original ar- 
ticles @J§]. Figure 1 shows a schematic version of the 
layered system. 




It has been argued from the theoretical point of 
view, |^-||] that lamellar and other microphase struc- 
tures can occur in mixtures of oppositely charged poly- 
electrolytes. Usually, such phases are associated with 
block copolymers where chemical bonding is an imped- 
iment to macrophase separation. In charged polymer 
mixtures, the entropy of mixing of the counterions for 
each of the homopolymer species provides a huge barrier 
to macrophase formation, allowing instead for the pos- 
sibility of microphases. In both cases, copolymers and 
polyelectrolytes, the separation is driven by the chemical 
repulsion between the species. 

An interesting aspect of the above outlined procedure 
for manufacturing the lamellar structure is that it would 
seem to work even at temperatures or for ranges of pa- 
rameters at which, in principle, a phase separation should 
not be observed. In the following we will analyze the 
conditions under which we can expect to find such sta- 
ble morphologies. Above the transition point, an RPA 
approach will show that a finite number of these layers 
should be stable near the surface. 

If the charged surface is immersed in a mixture of ho- 
mopolymers, it will drive the spontaneous creation of a 
few number of layers that will then coexist with the bulk 
of the mixture. This case is slightly different to the one 
resulting from the dip-and-wash mechanical layering pro- 
cedure in which the surface is only in contact with solu- 
tions of homopolymers. Both systems, however, should 
have similar stability conditions, the same relevant pa- 
rameters, and generate comparable length scales for the 
layers. We will carry out our calculation for the case of 
the surface immersed in the polymer mixture, and con- 
sider the charged surface clS cl el large flat plate. 



II. RPA FOR POLYELECTROLYTE MIXTURES 



FIG. 1. Scheme of the layered structure formed by 
a mixture of polyelectrolytes near a flat surface. The 
surface is charged positively, attracting a negatively 
charged polymer. The second layer consist of poly- 
mers positively charged, attracted to the first charged 
layer, and so forth. Counterions are not shown in the 
scheme. 



The fundamental fields for the description of the model 
system are the volume fraction of the homopolymers 
$,4 and Both homopolymers will be considered to 
have the same monomer number N . A fraction / of the 
monomers in each chain carry charges za and zb respec- 
tively and we will again consider only the most symmet- 
ric case z = za = —zb- We consider the polymers to 
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be dense with a negligible amount of solvent, and con- 
sider the effective volume fraction occupied by the rela- 
tive counterions very small. The local counterions densi- 
ties are ua and jib- Counterions for polymers A and B 
carry charges — za, and — zb- 

The RPA approximation for a mixture of polyelec- 
trolytes has been developed and studied by a number 
of groups H~||] ■ With a number of further simplifications 
(to be described later), a particularly useful form of the 
RPA free energy functional has been written by Nyrkova 
et al. To their basic free energy, we add a term cou- 
pling the charge density to an external potential (that 
will be generated by the charged surface). 
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In this expression, all energies are measured in units of 
ksT, and the Kuhn length is b. The Bjcrrum length 
is Bj = e 2 /(efcsT). The interaction parameter \ is m 
this case positive, implying repulsion between species. 
The parameter x is considered close to, but smaller than, 
the critical value required for the spontaneous creation of 
lamellar (or another) microstructure. The prefactor for 
the gradient terms is such that if we turn off the electro- 
static interactions, the correlation functions will match 
the first terms of the Debye function for each type of 
chain, namely, in Fourier space {^k'f-k} = C(l + + 
. . .), with the radius of gyration R 2 = b 2 N/6. Lastly, the 
density of charged particles p is defined by 



za(Ja^a ~ n A ) + z B (fB^B - Zb)- 



(2) 



Thus, we have essentially a regular solution model in 
which Nyrkova et al. J^] argue that there are no important 
corrections arising from fluctuations. Further, the intra- 
chain correlations are reduced to pure random walks, and 
only the first derivatives of the fields appear in the func- 
tional. In reference || the full Debye function has been 
used to represent the intra-chain correlations, but the 
effect on the results here presented is minor, as the cor- 
rections appear at wavelengths smaller than the resulting 
layer thickness. 

Since we are interested in the behavior of the mix- 
ture above the demixing temperature, we assume a base 
state with homogeneous distributions of homopolymers 
and counter-ions. The overall concentrations of polymers 
&A = &b = 1 — $0 are considered as given, and the 
bulk concentrations of counter- ions are = for each 
species. These overall concentrations should remain con- 
stant and the only important collective variables are the 
difference in the local concentrations between the species: 



* = $ A - $s - (2$o - 1) 
n = ru-n B -/(2$o-l)- 



(3) 
(4) 



Since the volume fraction of the counterions is small com- 
pared to the volume fraction of the polymers, we set 



<S>a + ®b = 1, 



(5) 



to be satisfied everywhere. 

The effective RPA free energy in terms of the n vari- 
ables is obtained by expanding all terms up to second 
order in the fields, and the result is 



F = F + ^ [ drdr ,(mr)-n(r)Kmr<)-n(r>)) 
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-Jdr^ 2 + J drV(r)ez(f^(r) -n(r)). (6) 

For our purposes, we need to Fourier transform all quan- 
tities (^(k) = J dr e _ik r 5'(r), and so forth) and rewrite 
the free energy as: 
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III. RESPONSE TO EXTERNAL POTENTIALS 
AND SURFACE EFFECTS 



For a given external potential V(r), we can now calcu- 
late the general response of the system. We shall work 
out the response when the external potential originates 
from a pair of parallel charged surfaces, and then, in more 
detail, the case of only one bounding surface. 



A. General effect of an external potential 

Minimization of the free energy, eq.|^, with respect to 
the polymer and counterion densities, gives a pair of cou- 
pled equations for ^ and h. To present the solutions in 
a more transparent manner, we will introduce a number 
of simplified and dimensionless variables, namely: 
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s = — [l- 2 X N<S>(1 - $)] 



t = 



4tt$(1 - $) Sj 
9 V' 



(8) 

(9) 

(10) 
(11) 



The variables s and t measure the strength of the inter- 
species and electrostatic interactions, while the variable 
x compares the wave- vectors against the the screening 
length of the counter- ions r =6(18/) -1 / 2 . 
Now the solution to the system reads 



¥(k) = JZi(k)V(k) 



[2a;$o(l-$o)](- e F(k)) 
D(k) 



(12) 



fl (k) =fl2 (k)V(k) = ^ 1 -^ + ^ e ^ (13) 

U(k) 

where the denominator is 

D(k) = x 2 + {t + s)x + t{l + s), (14) 

and we have introduced a shorthand notation for the 
overall responses R\ , and i?2 • 

B. Response of a mixture bounded by two surfaces 

We consider now a bulk of the mixture contained 
between two parallel surfaces. One of the surfaces is 
charged positively, with a surface charge density a. The 
more realistic condition for the second surface is to take it 
as uncharged. In fact, this second surface can be an imag- 
inary surface coinciding with the free surface of the mix- 
ture (in contact with, say, air). The amount of charge in 
the second surface determines the electrostatic potential 
difference between the two surfaces, which is the relevant 
driving force of the micro-segregation when we neglect 
non-linear terms. Thus, for simplicity, we consider in- 
stead the equivalent case in which one has two oppositely 
charged surfaces with equal absolute value of charge den- 
sity cr/2, separated by a distance L. In the following, we 
take the z axis to be normal to the surface, and set the 
origin in the left surface. 

From elementary electromagnetism, the potential cre- 
ated by the charged surfaces is, for values of the z coor- 
dinate between and L 



V(r) = V{z) 



,L 



z). 



(15) 



To use our Fourier transformed quantities, we extend the 
potential to the full z axis, by reflection at all planes with 
z = nL, with n an integer. This allows us to obtain a suit- 
able boundary behavior for the field and reduces the 
problem to consider a Fourier Series for the potential in 



the region (—L,L) (using the reflection condition). The 
expansion for the potential is V = V m e vm7TZ ^ L , with 
coefficients 



V m 



m odd 



2oL 

m even. 



(16) 



for The linear response to this potential is also expressed 
as a Fourier Series, and the coefficients for the the com- 
position ^ and the counterion density n are obtained 
immediately from the solution in eq. [l^ by evaluation at 
the wavenumbers k m , (k x — rmr/L, k y = 0, k z = 0): 

* m = R 1 (k m )V m , (17) 
n m = R 2 (k m )V m . (18) 



C. Response to the presence of one surface. 

The case of only one surface bounding the mixture is 
equivalent to the case considered in the previous sub- 
section when we let the separation between the surfaces 
grow and effectively become infinite. The proper limit to 
consider for the Fourier transform of the potential is 



V(k) - 2a(27T) 2 S{k x )S(k y )^. 



(19) 



It is convenient to introduce a new parameter propor- 
tional to the charge density, 



7 = Aear 2 ${l - $), 



(20) 



which amounts to measure the surface charge density as 
the number of charges in a square area of side r Q . 

In the final result, all quantities depend only on the z 
coordinate, and all relevant reference to the wave- vectors 
is only through x = k 2 /18f. We obtain that the trans- 
form of the response in this limit can be expressed as 



*(k) = (2ir) 2 6{k x )5(k y 
h(k) = (2ir) 2 6(k x )5(k 
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x 2 + (s + t)x + t(l + s) 
y> x 2 + {s + t)x + t(l + s) 



, (21) 
• (22) 



The denominator in these expressions is a real, even, 
quartic polynomial in |k|, and real quadratic in terms 
of x. 

In this mean field approximation, there is a transition 
to a micro-segregated phase when the roots of the poly- 
nomial become real. This occurs when, for a given value 
of t, the s variable becomes negative enough, namely, at 



-2t 1/2 + t, 



(23) 



and at that point, there are two pairs of double roots for 
the polynomial at 
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x r = ±{t 1/2 -t). (24) 

We are interested in the region above (temperature- 
wise), but close to the transition point. In such a region, 
the susceptibility to the influence of the surface is high, 
and yet, when the surface disappears there will not be 
spontaneous micro-segregation. We can measure the dis- 
tance from the transition point by the variable 

e = (s r - s) 1 ' 2 . (25) 

In terms of this parameterization the response has a 
rather simple expression. In real space we find that the 
field if? consist of decaying oscillations away from the 
z = plane, (i.e. the physical charged surface): 

mz) = —\ exp (-—z) cos f ^-^] . (26) 
e^r V r o J \ r J 



IV. LAYERING 

In the equilibrium compositional profile that we have 
obtained, eq. ^6], a layer is represented by a peak (pos- 
itive or negative) in the concentration variable. For a 
layer to be recognizable as such, the peak in the \& field 
has to be comparable in magnitude to 1 (in the case of 
symmetric mixtures) so that one of the components sat- 
urate a region in space thus forming a recognizable layer. 
Since the solution we have decays away from the surface, 
the number of layers that are formed is finite. It is clear 
that it is of interest to maximize the number of layers 
that are stable near the surface. We shall consider three 
ways, suggested by our results, in which this effect can 
be achieved. 

First, we note that that the field amplitude decays over 
a characteristic length T — ro/e. If we tune the charge 
at the surface so that the strength of the electric field 
is just enough to saturate the composition at the sur- 
face so that if? = 1 there, then layers will exist only up 
to a distance T away from the surface. The thickness 

of the layers is given by the wavelength of the oscilla- 
l li 

tions A = ro/27rx r , and therefore the number of layers 
formed near the surface is 

1 = ZttXj^ — ^ 
A e 

To acquire the largest number possible of layers in this 
way, we would like to make e as small as possible. Thus 
one has to be as close as possible to the transition point 
for the number of layers to suitably increase. An impor- 
tant issue arises here, in that while the layering system 
has been observed, the actual transition has not been 
observed, with the most simple explanation available be- 
ing that the parameter region at which the system will 



undergo a transition is out of the reach of physically real- 
izable systems, but then, even the near-transition region 
(e small), is as well hard to reach. 

A second idea for the creation of a good number of 
layers is over-saturation. Since we are interested in at- 
taining the *$> = 1 condition repeatedly, consider the case 
in which the source strength (the surface charge) is so 
big that, mathematically, if? will have to achieve values 
well beyond 1. To make sense of this solution, let us 
neglect non-linear effects and simply assume that when 
the mathematical solution dictates a value for "J larger 
than 1 the physical system simply saturates and takes 
the value of 1 in that region. Then, even if we go to 
distances well beyond T, the condition '5 = 1 will be sat- 
isfied many times away from the surface. To obtain m 
layers we need an amplitude e mA / r , so that 

m~log( 7 / 7 o), (28) 

where 70 is the effective surface charge density required 
to obtain saturation of one component at the surface, i.e. 
if! = 1. Schematically, we have the case shown in Figure 
2. While this is obviously a rather crude approximation, 
we should point out that, at the linear response level, the 
thickness of the layers does not depend on the strength 
of the surface charge density and we can expect that the 
layers formed might still have roughly such a character- 
istic size. A proper theoretical treatment of the problem 
should take into account the effects of strong segregation 
caused by the surface. 
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FIG. 2. The polymer composition profile 'J', away 
from a charged surface. When the surface charge is 
large enough, the amplitude as obtained in eq.Efjl ex- 
ceeds 1. Constructing profile for the composition by a 
simple truncation of this solution at the limiting values 
±1, we obtain important compositional fluctuations, 
even relatively far away from the surface. 

It is also important to note that with this kind of ar- 
gument we could attempt to reconcile the experimental 
observation of the layering process, with the fact that at 
the same time there is not know trace of the non-forced 
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microphase separation. This is, even when one operates 
far away from the critical point, a sufficiently large driv- 
ing force can show the traces of the micro-segregation. 
For simplicity, we have presented our results in a form 
valid only in the neighborhood of the transition. For val- 
ues of the parameters moderately distant from the tran- 
sition point case, one should use the full response, and 
look at the complex roots of the denominator D to ex- 
tract the wavelength, amplitude and decay length of the 
oscillations. 

Finally, we return to finite size effects. A third idea to 
enhance the number of stable layers is to place the sys- 
tem in a finite geometry, as was described in section III. 
This is roughly equivalent to the experimental situation 
in which the layered system is taken away from the bulk 
of the mixture and has an overall finite thickness. 

In the exploration of the limit when only one surface 
is present we have noticed that the oscillations of the ho- 
mopolymers concentration decay relatively quickly. This 
needs not to be the case when we have a finite distance 
between the surfaces. First, if the distance between them 
is roughly about twice the decay length L ~ 2r, we will 
have, in the non over-saturated case, twice the number of 
layers expected from our estimate eq. because of the 
association of the layers with each of the two surfaces. 
Further, the (possibly) constructive interference of the 
oscillations will make the amplitudes larger around the 
mid-point between the surfaces, creating a stronger seg- 
regation pattern. Secondly, we can try to match a mul- 
tiple of the preferred wavelength mX with the distance 
between surfaces. Clearly, this resonance maximizes the 
segregation of the species. In Figure 3 we show the re- 
sponse constructed from our result eq. [l7], for a distance 
between surfaces equal to 5.5 times the wavelength A (9 
layers). 




FIG. 3. The polymer compositions profile 

^ for a finite slab of polymer mixture of width 
L — Hiiro/xl^ 2 — 38ro. The values of the other basic 
parameters are t = 0.5, s — —.91. 



V. DISCUSSION 

Our description of the response to the charged surface 
has been presented in terms of a set of reduced variables, 
ro,s, t. But it is necessary to comment on the values 
these quantities can achieve in actual systems. The sim- 
plest quantity is t. This is just a multiple of the Bjerrum 
length and we have it will always oscillate around a value 
of 1, taking values of, in some cases, 1/2 for stiff poly- 
mers. 

Next we must consider the values of ro- This is per- 
haps the more important consideration in the problem. 
We note that the root of the determinant x, in the vicin- 
ity of the critical point x r , is roughly of the order of t, 
and from our previous paragraph, also of the order of x r . 
Thus, in this region, s will take values of order 1. The 
oscillation frequency is then always of the same scale as 
the screening length. Typically we will be interested in 
creating layers of thickness comparable, at the very least, 
to the radius of gyration of a chain, R 2 ~ N, which to- 
gether with rg ~ 1//, leads to the condition 

fN ~ 1. (29) 

For this condition to be satisfied it is required that there 
should be only about 1 free counterion per chain, or at 
least a very small number of them. The reason while 
this condition might still be met in practice is the possi- 
bility that many of the counterions are tightly bound to 
the chains in the melt, and / represents only an effective 
number of free counterions. 

The reported findings of the Decher group , seem 
consistent nevertheless with this picture. While they ob- 
tain mono-layers, the thickness of the layers is relatively 
small, about 100^4°, and the observed layering occurs also 
under weakly charged conditions. 

In the neighborhood of the critical point of the mixture 
the parameter s is also of order — 1, which implies that 

X*f~fN~l, (30) 

which simply implies that the segregation forces should 
be strong enough to compensate for the loss of transla- 
tional entropy of the effectively free ions. This is not a 
stringent requirement and many known systems can sat- 
isfy the condition. 

Consider, finally, the charge density for the surface. 
Even when the systems is not too close to the microphase 
segregation point, and e is not too small, the amplitude 
of the oscillation can be set to 1 or be over-saturated by 
increasing 7. If we put q charges in a surface area of 
size R 2 ~ N, we will have 7 ~ qj fN ~ q. It is then 
clearly feasible to charge the surface with a large number 
of charges q in the area occupied by one chain. The ap- 
proach of obtaining a large number of layers using a very 
strong driving force seems then easily realizable. 
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VI. CONCLUSIONS 
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We have explored the possibility of creating stable mi- 
crostructures in polyelectrolytes by means of a strong 
external influence. While our qualitative results are con- 
sistent with experimental work, a more precise compari- 
son with experiments cannot be done directly. We have 
used a number of effective parameters some elucidation 
is required to relate them to the bare parameters of the 
system. For example, it is necessary to develop a relation 
between the original number of charges in every chain 
with the effective number of free counterions, and also 
to the effective value of the \ parameter, which should 
also be affected as the number of condensed counterions 
changes. These are subjects of forthcoming research. 

The only property that we have specified for the sur- 
face is its charge density, but is clear that it can be made 
of different materials, and in some cases the surface can 
already have a previous polymeric coating. The layering 
scenario can then occur, for example, in the neighbor- 
hood of surfaces of charged polymer brushes jfi| . 

Another aspect of the problem that should be mention 
is the possibility of forcing different structures other than 
lamellar-like layers. When the volume fraction for one of 
the components is small, one could expect that, as is 
the case in block copolymers, the micro-segregation will 
create other structures, such as spherical and cylindrical 
domains. Of course, there would also be an interplay be- 
tween the geometry of the external charged surface and 
the preferred geometry of the segregated system. 
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